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We study single-particle excitations in the Mott insulator phase of a Bose gas in an optical lattice. 
The characteristic feature of the single-particle spectrum in the Mott insulator phase is the existence 
of an energy gap between the particle and hole excitations. We show that this energy gap can be 
directly probed by an output coupling experiment. We apply the general expression for the output 
current derived by Luxat and Griffin, which is given in terms of the single-particle Green's functions 
of a trapped Bose gas, to the Mott insulator phase using the Bose-Hubbard model. We find that 
the energy spectrum of the momentum-resolved output current exhibits two characteristic peaks 
corresponding to the particle and hole excitations, and thus it can be used to detect the transition 
^ , point from the Mott insulator to superfluid phase where the energy gap disappears. 
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I. INTRODUCTION 



'"^ , Ultracold atomic gases in optical lattices provide a new framework for the experimental study of many-body 
■ quantum systems by making use of their remarkable degree of experimental control of the relevant parameters. In 
O I particular a recent pioneering experiment by Grcincr et al. P| observed a phase transition from the superfluid to Mott 
^ . insulator phase in a ultracold gas of *^Rb atoms trapped in a three-dimensional simple-cubic optical potential. 

In the experiment of Ref. Q , the characteristic feature of the superfluid phase appeared in the interference pattern in 
CsJ , the absorption image of the atomic cloud after ballistic expansion, which is due to the coherence of the Bose condensed 
^ ■ atoms. The phase transition to the Mott insulator phase was signaled by the disappearance of the interference pattern. 
, However, it is not obvious how one canprecisely identify the transition point from the interference visibility. Moreover, 
' recent papers by Roth and Burnett [3, pi argue that the interference pattern contains only the information on the 

^ ] quasi-momentum distribution of the lattice system, but has no direct information on the superfluid fraction. Therefore 

. the disappearance of the interference pattern does not provide direct proof of the Mott transition |^ Q . On the 
' other hand, in the Mott insulator phase, their characteristic feature appears in its single-particle excitation spectra, 
, which have a gap between particle and hole excitations. The experiment by Greiner et al. also provided the 
' spectroscopic study of this energy gap by applying a potential gradient. 

, In this paper, we propose to use an output coupling current as a more direct probe for detecting the Mott gap of a 
^ ■ Bose gas in an optical lattice. Recent papers ^ JJj discussed the output coupling current from the trapped atoms 
\ at finite temperature, and pointed out that the analysis of the output current serves as a probe of many-body states 
of a trapped gas. In particular, Luxat and Griffin Q developed a general formalism of the output coupling current 
Q I by using the linear response theory. Their formalism is useful because one only has to concentrate on calculating 
^ I single-particle correlation functions of the system of interest. In this paper, we calculate the output coupling current 
from a Bose gas in an optical lattice in the Mott insulator phase, and explicitly show how the characteristic feature 
of single-particle excitations appears in the energy spectrum of the output-coupling current. 

To describe a Bose gas in an optical lattice, we use a Bose-Hubbard model [a, Si- This model has been shown to 
J-j \ exhibit the phase transition from the superfluid to the Mott insulator at zero-temperature, and its phase diagram 
. - 1 has been investigated by many authors both analytically and numerically |3, ITol IllL IT^ IT^ IT3 . ITsI l . The mean-field 
theorypredicts the transition point, at which the order parameter for the superfluid states vanishes, U / zt = 5.83 = 
Uc/zt [sllllllTsll. where U is an on-site interaction, i is a hopping parameter, and z is the number of the nearest sites 
in a Bose-Hubbard model. More recently, finite-temperature behaviors of a Bose gas in an optical lattice have been 
studied in Refs. In Sec.|nl we briefly discuss the mean-field phase diagram of a Bose-Hubbard model at finite 

temperatures using the Landau theory of the phase transition. A detailed microscopic derivation of the Landau free 
energy of a Bose-Hubbard model is presented in Appendi xjXl 

In Sec. mil we use the formalism of Luxat and Griffin [3 to express the output coupling current from an optical 
lattice in terms of the single-particle correlation functions of the Bose-Hubbard model. In Sec. lIVI we then calculate the 
sin gle- particle Green's function of the Bose-Hubbard model by using the standard basis operator formalism developed 
in [13]. This formalism is useful for calculating excitations in the strong-coupling ([/ > t) regime near the Mott 
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transition. We solve the equation of motion for the single-particle Green's function in the Mott insulator phase by the 
random-phase decoupling approximation. The poles of this Green's function provide dispersion relations of particle 
and hole excitations. These dispersion relations exhibit an energy gap, which vanishes at a critical point and grows 
with increasing on-site interaction as well as with increasing temperature. We will show that this behavior of the 
energy gap can be directly probed by the output coupling current. 

In Sec. we give an analytical expression for the output coupling current from an optical lattice in terms of the 
spcctrul functions of the Bose-Hubbard model calculated in Sec. lIVI We then explicitly calculate the output current as 
a function of the energy transfer from probe laser field, and show that a momentum-resolved output coupling current 
exhibits two characteristic peaks corresponding to the particle and hole excitations, allowing for direct observation 
of the Mott gap. Thus, one can use a momentum-resolved output-coupling current as a direct probe of the phase 
transition from a Mott insulator phase to superfluid phase. 



II. MEAN-FIELD PHASE DIAGRAM OF A BOSE-HUBBARD MODEL 



In this section, we discuss the superfluid-Mott insulator transition of a Bose gas trapped in an optical lattice 
potential. Our system is described by the Hamiltonian 



= J dr 4ir) 
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^— V2 + Kp(r)-f yt(r)-/i 
dr '0|(r)?/>J(r)'0t(r)V't(r), 



(1) 



where tpti^) is the Bose field operator for a trapped atomic state. Here Vt{r) is the magnetic trap potential and 
K)p(r) = Vb X]j=i sin^(27ra:j/A) is the simple-cubic optical lattice potential with the lattice constant a = A/2. The 
interatomic interaction is approximated by the short-range potential with the s-wave scattering length Og. We assume 
that the band gap energy between the first and the second excitation band is large compared to the temperature, and 
thus only the first band is thermally occupied. The depth of the optical lattice wells is assumed large enough to make 
the atomic wave functions well localized on the individual lattice sites. Under these assumptions, one can expand the 
field operator V't(r) of the trapped atoms in terms of the Wannier function w{r) 



(2) 



and rewrite Ht into the Bose-Hubbard model 



<ij> i i 

where b] and bi are the Bose creation and annihilation operators of atoms at the lattice site z, respectively. The 



parameters in Eq.Q are defined by 

t = J dr w*{r — Ti) 

dr \w{r)\\ 



wir-rj), 
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= fi- J dr Vt(r)|w(r - r, 



(4) 
(5) 

(6) 



where t is the hopping parameter for the adjacent sites i and j, U is the on-site repulsion, and fXi is the site dependent 
chemical potential, given by the sum of the chemical potential /i and the energy offset e,; of each lattice site due to 
the slowly varying magnetic trap potential Vt- In Eq.(j2J), the summation in the first term on the right-hand side is 
restricted to nearest neighbor lattice sites. 

In this paper, we ignore the confining magnetic trap potential Vt, setting = /i in the Bose-Hubbard model 
As mentioned in Sec. P the properties of the Bose-Hubbard model have been investigated extensively by many 
authors H, 0, 0, il5| . The effects of the confining trap potential to the equilibrium properties of the 
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FIG. 1: The phase diagram of the Bose-Hubbard Hamiltonian. The solid line is the phase boundary for 13 zt 
dashed line is for /3zt = 2.0. The Mott phase has the filling factor (average particle number per site) n = 1. 
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Bose-Hubbard model have been studied in Refs. [20ll2ll|. Here we discuss the superfluid-Mott phase transition in this 
system using the Landau free energy, which is the function of the order parameter in the superfluid phase defined 
by 4" = (bi). In Appendix ^ we give a detailed microscopic derivation of the Landau free energy by means of the 
inversion method developed by Fukuda et aL[l9l|. The resulting Landau free energy per site is given as 

r{-9,-^*)^r'-°^ + A\^\^ + B\^\\ (7) 

where the explicit expressions for F^''^ and the coefficients A and B are given in terms of the parameters U , t, and 
/i in Appendix 1X1 From the free energy in Eq.©, one finds = for ^ > and 'J 7^ for A < 0. Thus the 
phase boundary between the non-superfluid phase and the superfluid phase with ^' 7^ is found by setting ^ = 0. 
We note that our Landau theory at T = reproduces the mean-field results in Refs. ITllll5l|. Figure ^ shows the 
phase diagram for the temperatures I3zt = 10.0 and fizt ~ 2.0. The similar phase diagram has been also obtained in 
Refs. HQ 

Fi gure U shows that the region of the non-superfluid phase grows as the temperature increases. As discussed in 
Ref. [13, this does not mean that the region of the Mott insulator phase grows because the Mott insulator phase is 
meaningful only at zero-temperature. On the other hand, one can call the low-temperature non-superfluid phase at 
(izt = 10.0 a Mott insulator phase for practical purpose , since the compressibility is very close to zero. We also note 
that the phase boundary at (3zt = 10.0 is almost identical to that of T = 0. 

In Sec. we will investigate excitations of a Bose gas in the optical lattice, modeled by the Bose-Hubbard 
Hamoltonian in Eq.J^j), near the phase boundary with the flxed value of the chemical potential fi/ zt ~ 2.41 approaching 
from the Mott (or normal) phase by using an output coupling current. This is the main purpose of this paper. 



III. OUTPUT COUPLING CURRENT EXPRESSED BY CORRELATION FUNCTIONS 

In the output coupling experiments, one uses a stimulated Raman transition to change the internal states of atoms 
from a magnetically trapped state to an untrapped state, and then extracts atoms from a trap to a free space. This 
technique was originally introduced for the atom laser device in order to extract a coherent atomic beam from a 
Bose-Einstein condensate in a trap potential As pointed out in Refs. 0,0,01 this technique can also be used 
to probe many-body states in a trapped Bose gas. In the present paper, we consider the experimental situation 
of Ref. Q, where the atoms are initially trapped by using the combined magnetic and optical lattice potentials, as 
described by Ht in Eq.([3l. Although we ignore the effect of conflning potential in the actual calculations, we still 
consider the situation where the gas is confined in the finite region due to the magnetic trap potential. By using 
the Raman transition, one can change the internal state of atoms to magnetically untrapped state associated with 
momentum kick, and push atoms out of the confining region. In general, atoms in the output state will feel an optical 
lattice potential produced by the laser field. Following the notation of Ref. 0, for convenience we call the initial 
and output atomic states "trapped" and "free" states respectively, although atoms in both states feel optical lattice 
potentials. In this section, we express the output coupling current from a Bose gas in the optical lattice in terms of 
single-particle correlation functions of the Bose-Hubbard model, following the formalism of Luxat and Griffin . 
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The total Hamiltonian of our system is given by 



n = Ht+Uf + Vit) 
= H° + V{t), 



(8) 
(9) 



where Tit describes the motion of Bose atoms in a combined magnetic and optical lattice potential and is reduced to 
the Bose- Hubbard model in Eq. ^ . The second term of Eq. (jSJ describes the motion of atoms in a free state : 



Hi = J dri'lir) 

dr V't (r) [hiir) - ^f] 



^f(r) 



(10) 



where V'f (r) is the Bose field operator, fif is the chemical potential, and V^p{r) is an optical potential for atoms in this 
state 123 • We assume that the density of a Bose gas in the free state is so low that the interaction between atoms in 
the gas can be neglected. The third term of Eq.® describes the time-dependent probe laser field 0: 



V{t) = 7 / rfr 



„i(qT-wt) / t 



^T(r)V^t(r)+H.c 



(11) 



where huj is the energy transfer and hq is the momentum transfer from the laser field. We assume that the trapped 
system and untrapped system do not interact with one another without V(t). The time dependent external field 
V(t) couples two systems, acting like the tunneling Hamiltonian. This field annihilates an atom in one system, while 
creates an atom in the other system and vice versa. In this process, the laser field causes transition between two levels 
and then one can observe the extracted atoms. 

We define the output coupling current from a gas in the trapped state to a gas in the free state 



dt ' 



(12) 



where Nt — J dripl{r)ipt{r) and (• • •)t indicates a nonequilibrium statistical average. The negative sign in this 
expression implies that the current from a gas in the trapped state to a gas in the free state is positive. By using 
the linear response theory regarding the time dependent external field V(t) as a perturbation, tunneling current H12|) 
becomes 



61 = T^Re y dr y dr' e^'^ ^"-'-') ^ ^ 



C^l^, (r'r, -^')C^,^t (r'r, + ^ - A^/?.) 



with the correlation functions defined by 

C^t^(rr',0 = (^\rt)i:{r\t^O)) 
= Tr{peqV^t(rt)^(i.')} 

= iTr{e-^«°e'«°*/V^(r)e-*^''*/V(r')} 
'z = Tr e-^'^" 



(13) 



(14) 



where 



-f3H" 



is the initial equilibrium density operator in the absence of the external field, ^j{t) 



-iHot 



/^Tpe ^''''■ot/h i-j^g Held operator in the interaction picture, and A/i = jif — fi. The Fourier transform of 



the correlation functions is defined by 



C^tv,(rr',i) = r *C^tv,(rr',a;). 



(15) 
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Using the correlation functions, we define a single-particle spectral function as 

/oc 
dt e-*([V'(rt),^t(r')]) 
-OO 

= C^^,t(rr',cj) - C^t^(r'r,-cj). (16) 

One can write the correlation functions C^t^ and C^^t in terms of this spectral function as 0,123 

C^t^(r'r,-tj) = fihu;)A^,f^{r'r,uj) (17) 
C^^t(rr',c^) = [l + /(M]^^tv,(r'r,t^), (18) 



where f{Tiuj) is the Bosc distribution function 



/(^-) = ^^r^i (19) 



The spectral functions for the non-interacting atoms in an optical lattice can be expressed in terms of the eigen- 
function (or Bloch function) Up(r)e*P '" and eigenvalue ef(p) of the single-particle Hamiltonian hf 

A^t^^(r'r,^) = 2^^Up(r)u;(r')e-^P-('-'-'-)j(^- [ef(p) -A*f]A), (20) 
p 

where p is the quasi-momcntum of the Bloch state. On the other hand, the correlation functions for the trapped 
atoms are rewritten by using the Wannicr functions: 

C^t^^(r'r,t) = {4{v't)U^)) 

= Y.^*{v' ~Vk)w{v-vi){hl{t)h). (21) 

In order to calculate the correlation functions for the trapped atoms, we must calculate the time-correlation functions 
such as {b\{t)bi). 

For calculating the output-coupling current (|13|l . we need the following four correlation functions: 



/OO 
dt e-'^''{bi{t)bi), 
-OO 



(22) 



C, ,t(rr',w') 

/OO 
dt e'^''{bk{t)b\), 

= 27r^ e-*P-<'-'"') [1 + /(W + hLo- A/^)] 
p 

X up(r')w;(r)5(w' +UJ- [ef (p) - y]/K), 
p 

X Up{r')Up{r)S{uj' + uj - [ef (p) - ^i]/h). 



(23) 



(24) 



(25) 



The purpose of the next section is to calculate the time-correlation functions (6[.(t)&;) and {bk{t)bj) of the Bose- 
Hubbard model that appear in the expressions for C'^t^^ and C'^^.^t given in Eas. H22|l and H23|) . 
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IV. CORRELATION FUNCTIONS OF THE BOSE-HUBBARD MODEL 
A. Standard-basis operator method for Bose-Hubbard model 

In this section, we apply the standard-basis operator method dev elop ed in Ref. to calculate the single-particle 
correlation functions of the Bose-Hubbard model. We note that Ref. ^l^ also used this formalism to discuss excitations 
in the Bose-Hubbard model. Here we give detailed calculations of the correlation functions of finite temperatures. 
The standard-basis operator is defined as 

V^^,^\ia){ia'\, (26) 

where a complete set of functions for i = 1, 2, TVgite and a = 1, 2, consists of the state vectors of atoms 

at the site i in state a. The set of states chosen is arbitrary. Any single-site operator O then can be expressed as 

O = ^^\ia){ia\0\ia'){ia'\ 

i a., a.' 

= ^^(^a|0|za')iLa'- (27) 

i a, a' 

Using the standard-basis operator, one can rewrite the Bose-Hubbard Hamiltonian in Eq.® as 

= - X! X!X!^"a',/3/3'-^aa'^^/3'+X!X!^"«'^aaM (28) 

<ij> a, a' /3,/3' i a, a' 

where 

Ki-,f30' ^ U,{^a\bl\^a'){JP\b,\JP') 

+ Uj{ta\b,\za'){jp\b]\jP'), (29) 

^aa' = ^{^c^\n'\^c.') - + ^) (30) 

We define the two-time retarded and advanced Green's functions in the coordinate representation in terms of the 
standard-basis operator by 

^ ~z0(t-t')([^L'(O,i|,/3'(i')]>, (31) 
GiV(^-i') = ((iL'(i)|iV(^')»A 

EE ^9{t'~t){[V^^,{t),Ll^,{t')]). (32) 

where 9{t) is the Heviside step function and [ , ] denotes the usual commutation relation. One can interpret the 
retarded Green's function as the probability amplitude that at time t' the jth site atom transits from the state to 
the state /? then at time t the ith atom transits from the state a' to the state a. The Fourier transforms of these 
Green's functions are defined by 

dt e-(*-*')Gfifj3|(t-t')- (33) 

As discussed in Ref. [S^, one can show that the retarded (advanced) Green's function G^{uj) {G^{uj)) can be ana- 
lytically continued into the region of complex lo in the upper (lower) half-plane. If we make a cut along the real axis, 
the function 

''-''''^^-\Gtj^,,i.) for lm.<0. ^''^ 



can be regarded as a single analytic function, consisting of two branches, one of which is defined in the upper and the 
other in the lower half-plane of complex values of lj. 
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The above Green's functions are related to the time-correlation function through [2S 



X /(;ia')e-''"(*-*') . 



(35) 



where S is an infinitesimally small positive quantity. The relation (|35|) is essential for calculating the time-correlation 
functions in and if^ . 

Differentiating the Green's function (|33|l with respect to t and Fourier transforming into uj, we obtain the equation 
of motion 



(36) 



This equation is not closed for G^^, ^^/(w) since it involves the higher-order Green's function, i.e., three-operator 
Green's function arising from (([L^^,,7it]|i^^/))w Therefore we need to decouple these terms in terms of the two- 
operator Green's functions. For this purpose, we first rewrite Ea. H36(l by using a commutation relation for 



The equation of motion then becomes 

I 11 vv' y 



(37) 



(38) 



In the random-phase approximation (RPA), the three-operator Green's functions are decoupled by two-operator 
Green's function as llSll 



- Sc.,{LUG'j,,^pf,,iu;) + ^..'(iL)GLV,^^,(^). 
Substituting Eqs.JSHIl and (gOJ into Eq.(|nHl), one obtains 

^^G^^„,^^^,(u;) 



(39) 
(40) 



E 



{y^', -EE Din',^,AG^^^^,,,{co) 

I v ^ 



(41) 
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where represents the probabiUty that the atom at the site i is in the state a, which is defined by 

K ^ {L'U ■ (42) 

Hereafter we consider a homogeneous system ignoring tlie trap potential . i.e., fii = ji. In this case, Ea. (|41fl can be 
transformed into the k-space: 

?iwGQ,a',;3/3'(k, 
= Sa/3'Si3a'h{Da - Da') 



E 



where 



-'Vsitc , 



kelstBZ 



with iVsito being the nmnber of the lattice sites and 

rpk \ ^ rjilm ?k-(ri— r^) 



(43) 



(44) 



(45) 



The summation of the wave vector k is restricted to the first Brillouin zone. In order to solve Ea. (|43|l . we must specify 
appropriate basis {Iza}}. 



B. Green's functions 



The time-correlation functions in H22|) and (|23|l are expressed in terms of the standard-basis operators as 

aa' /3f3' 

= EE(*"i^^N"')0'/3i^^ij/3')(iLo'Wi'^./3')- 

aa' /3/3' 



Using the relation we find 



= Y.Y.^^a\b\\^a'){J(3\b^J|3') 



aa' 13/3' 

duj 



G^P'.aa'(^-^S)-G^^,^^^,(^ + ^S) 



m)b]) 



(46) 



(47) 



(48) 
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aa' 0(3' 

J-oo 27ri 



(49) 



For solving Ea. (|43|) in the Mott insulator phase, we take the Fock states {\n)} as the complete set. In order to obtain 
a closed set of equations for the Green's functions, we restrict the Hilbert space to {|n — 1), \n), \n + 1)}, where n is 
the number of boson per lattice site. We note that these states are used in the second-order perturbation calculations 
in terms of the hopping term at T = With this restricted Hilbert space, Eqs. H48|l and (|49|l become 

= ^^(z«|6l|za')0-/3|&,|j/?') 

aa' 013' 



X [Gppi^aa' (k, LU -iS) - Gpp'^aa' (k, LU + iS)] 

y r — e»{-*-k-(r.-r.)} f^r^) 
^'^^^ kGlstBZ ■'^'^ 

X [Gi (k, w - - Gi (k, uj + i5)] , 



{h{t)h]) 

= ^5](z«|&.|z«')(j/3|&]|j/3') 



aa' PP' 



(50) 



X [Gpp'^aa' (k, LO -iS) - Gfjp'^aa' (k, UJ + 2(5)] 
^^^^ kGlstBZ 

X {G2(k,u;-i5)-G2(k,c^ + 2,5)}, (51) 

where we have introduced the Green's functions Gi and G2 defined in terms of the Bose operators 

Gl^"(t-t') ^ {{b\{t)\b,{t'))), (52) 
G'i{t-t') ^ {{hm]{t'))). (53) 

and their Fourier transforms 



Gf(c.) = -^ ^ e-*''-(---'-^)GKk,u.), l = \,2. (54) 

-'Vsitc 



kGlstBZ 

These are related to the Green's functions defined in terms of the standard basis operators 

Gi(k,c.) = ^^(*«|6||*a')(j-/3|6,b-/3')G^/3',aa'(k,o;), (55) 

aa' PP' 

G2(k,c.) = Y.Y.^ia\h\ia'){jP\h]\]P')Gpp',^^,{\^,oj). (56) 

aa' pp' 

Recalling that we restrict our Hilbert space to the Fock space {|n — 1), \n + 1)}, these are expressed as 

Gi(k,w) = (n + l)G(„)(„+i)^(„+i)(„)(k,a;) 
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+ \/n{n + l)G(„)(„+i)^(„)(„_i)(k, 
+ \Jn{n + l)G(„^i)(„)_(„+i)(„)(k, 



(n)(«-l) 



(57) 



G2(k,w) = (?^ + l)G(„+i)(„)_(„)(„+i)(k, w) 

+ V^n(nTT)G(„)(„_i)_(„)(„+i)(k,tj) 

+ •v/n(n + l)G(„+i)(„)^(„_i)(„)(k,tj) 

+ "G(„)(„_i)^(„_i)(„)(k, cj), 



(58) 



where we omitted the site index of n as we consider a homogeneous system. Introducing the spectral functions for 
Gi and Gi 



Ai(k,w) = [Gi(k,cj-i(5)-Gi(k,w + i<5)], 
A2(k,w) = -i[G2(k,c^-i(5)-G2(k,w + i<5)], 

we obtain the correlation functions as 



(59) 
(60) 



kelstBZ 



i{wt-k-(ri-rj)} 



iVsit 



kGlstBZ 



dijj 



f{hio)A,{k,u;), 



(61) 



I 2^ 



,i{wt-k-(ri-rj)} 



/(nc^)A2(k,L^). 



(62) 

The explicit forms of the Green's functions H57|l and H58|) are obtained by solving the equation of motion (|43|l . The 
components for Gi(k, w) and G2(k, w) in Eas. (|57() and (|58|l are given in Appendix IbI The poles of these Green's 
functions determine the dispersion relations of the particle and hole excitations. These are given by 

U , 



ep,h(k) = —(271 - 1) - A* 



{n + l)P„+i,„ nPn,n-i 



6(k) ± -A(k), 



(63) 



where we denoted Dn 



Pajj and 



A(k) 



C/2 + 2U {{n + l)P„+i,„ - nPn,n-i} e{k) 

+ {2n{n + l)P„ + l.nPn,n-l + {n + lfPn+i,n 

-n^Pl,,_,}e{kf'" 



ep(k) - eh(k), 

3 

e(k) = 2t cos(fcjQ), 



(64) 
(65) 



with the lattice constant a = A/2. The dispersion relations in Ea. (|63|l depend on temperature through the function 
P. To obtain an explicit expression for the temperature-dependent function P, we must determine the occupation 
probability in a self-consistent manner. To do this, we follow the procedure described in Sec. VI of Ref. J^]. This 
procedure is rather involved, but the final results in the low-temperature limit _D„ ^ Dn±i are simply given by (see 
also Eq.(9.14) of Ref. [H) 



site 



J2 [/(6p(k))+/(|6h(k)|)] 



(66) 



kelstBZ 
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FIG. 2: The dispersion relations of the particle excitation ep(k) (upper branch) and the hole excitation eh(k) (lower branch) for 
k = (fc,0,0). Panel (a) shows the U dependence at the fixed temperature Pzt = 10.0, while panel (b) shows the temperature 
dependence at U = Uc = 5.83zt. 



— /(^p(k)), (67) 

^^^'^ kelstBZ 

^"-1 = 7J- /(kh(k)l), (68) 



^"+1 ^ M. 

kelstBZ 



kelstBZ 

where /(|eh(k)|) = -1 - /(-eh(k)), and 

Pn,n+i = 1 + E [l-2/(ep(k)) + /(eh(k))], (69) 

^'^^ kelstBZ 

Pn,n-i = 1 + E [2-/(ep(k)) + 2/(eh(k))]. (70) 

^^""'^ kelstBZ 

The right-hand sides of Eqs. (|69|l and (|70|l still depend on P through ep(k) and eh(k), and thus we must solve these 
equations along with (|63() self-consistently. In the very low temperature regime of interest, we can solve these equations 
by iterating (|69|) and H7U|) with H63|l . Hereafter we use the first-order results of this iterative approximation. 

Using Eas. (jBl|) - (jB4p in Eq.|S7|), we obtain the Green's function Gi and the spectral function Ai in the following 
explicit forms : 

r rir - ^^^^ I ^-^i^) ^7^^ 

ijj~e-p\K)/Ti w — eh(k)/n, 
Ai(k,tj) 27rZ(k)5(w-ep(k)/n) 

+ 2^(a-Z(k))<5(c^-e^(k)/n), (72) 

where the wave-function renormalization factor is defined as 

^(k) = ^i-y|(n + l)P„,„+i[^-[/(n-l)-?7P„.„_ie(k)] 
- nPn,n-i [fJ. - Un +{n + l)P„^„+ie(k)] 



+ 2n(n + l)P„,„_iP„,„+ie(k) + ep(k)|, (73) 

with 

a= {n + l)P„,„+i - nPn,n-i- (74) 
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One can also obtain analogous expressions for G2 and ^2- In order for Ai given in Ea. (|72() to satisfy the sum rule 
for the spectral function, a must be 1. At T = 0, one can show that P„,n+i = P-n.n-i — 1 ^^nd thus a = 1. However, 
our numerical evaluation of Ea. H74() shows that a slightly deviates from 1 as we increase the temperature. This is 
because we restrict our Hilbert space to {\n — 1), \n), \n + 1)}. In our calculations shown below, we only consider the 
low temperature region where a « 1 is satisfied (within 1%). 

At T = 0, the Green's function and the spectral function reduce to 

cj — ep(k)/7i uj — t\i[k.)/Ti 
Ai(k,w) =27rZ(k)(5(w-ep(k)/?i) 

+27r(l-Z(k)),5(t^-e„(k)/n), (76) 

where the wave-function renormalization factor at T = is 

and the dispersion relations for the particle and hole excitation at T = are given by 

6p,h(k) = |(2n-l)-^-^ 
1 



± 2 - + 2)C/e(k) + e(k)2. (78) 



These results for T = agree with results in Refs. [la, 113 • 

It is now useful to note the relation between the excitation energies and the phase boundary discussed in Sec. ^1 
The coefficient A in the Landau free energy in Eq.© is related to the static susceptibility to the symmetry-breaking 
external field (sec Appendix . and thus one has the relation 

^ (79) 



Gi(O,0) 



From the expression for Z(k) in Eq. we see that at the transition point A vanishes as ^ cx Aq, where Aq = 

A(k = 0) is the energy gap representing the difference between the particle and hole excitations at k = 0. Thus, 
Aq — signals the phase boundary between the superfluid and Mott (or normal) phase. In this paper, we have 
obtained A and Gi(0, 0) independently using the different approaches. At T = 0, it is straightforward to show that 
Ea. (|79|l holds by comparing Eq.jTSJ and Eq. (jA34l) . At finit T, we numerically confirmed that the transition point 
determined from A = Q agrees with the vanishing point of Aq. 

In Fig[3 we plot the dispersion relations given in Ea. (|63|l . With the fixed temperature (3zt = 10.0, the energy gap 
vanishes ai U — Uc — 5.832t and increases as U / zt increases (Figl^Ja)). In Figl^Jb) for the fixed U ~ Uc, we find the 
growth of the energy gap with the increasing temperature. These behavior of the energy gap is consistent with the 
phase diagram in Fig^ As discussed in Sec^ although there is no Mott insulator phase at finite temperature in a 
strict sense, we may call the energy gap Aq shown in Figl^a) for (3zt = 10.0 the Mott gap in a practical sense, since 
the compressibility is very low jig . 

V. OUTPUT COUPLING CURRENT 

A. General expressions for the output coupling current 

We can now express the output coupling current given in Eq. p3|l using the spectral functions derived in the previous 
section. The correlation functions (|22|l and (|23(l are written as 



J- — ^ w*{r' — ri)w{r — r^) 

'site . . 

^ e''^-("---^^)/(W)Ai(k,c.'), (80) 



X 

kGlstBZ 
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C, ,t(rr',w') 



*sitc . . 

J2 e'^<'''-'''^f{-hLu')A2{k,~Lu'). (81) 



X 

kelstBZ 



Using Eas. (|Sn|l and (|HTT) along with Eas. lf^ and in Eg. lIT^ and introducing the function ^^'(q, p) as 

0(q,p)^ / drw{r)u;{r)e^^^-P>'^, (82) 



we obtain the following expression for the output-coupling current: 

site kelstBZ p 

X 1^ (k + p - q) [1 + /(ef(p) - W)] /(ef (P) - [h^ + M])^i(k, - [cpn - ^^]/h) 

- (5 (p - k - q) /(ef (p) - /if) [1 + /(ef (p) - [hco + ^2(k, -c^ + [et(p) - /i]/?i)| . 

(83) 

In Ea. H83|) . the 5-functions represents the conservation of the total momentum. The first term in Ea. (|83|) describes 
a current of atoms going out from the combined trap, while the second term describes a current of atoms going into 
the combined trap. 

When we assume that /(cf (p) — Mf) ^ 1, which results in negligible tunneling from untrapped atomic gas back to 
the combined trap region, the above expression reduces to 

^^ = §^ E Ei^(q'P)i^ 

site kelstBZ p 

X ^ (k + p - q) /(ef (p) - [n^ + M])Ai(k, UJ - [ef (p) - fi]/h) (84) 

We have thus arrived at the general form of the output-coupling current from an optical lattice in the non-supcrfluid 
(i.e., Mott insulator or normal) phase in terms of the spectral function. 

B. Numerical evaluation 

We now calculate the output current given by Eq. (|84|l numerically. The sums over p and k in Eq. (|84|) can be now 
replaced with the integrals through V J rfp/(27r)^ and X^keistBZ ~^ ^ I-i"/a '^k/(27r)^, where V ~ NsUca^ is 

the volume of the system and a is the lattice constant. Then the output coupling current becomes 



/ 7^ dk|</,(q,p)p<5(k + p-q) 

Z(k)<5(ef (p)/;i - ep(k)/a -[u; + /i/;i])/(ep(k)) 



[a - Zik))S{ef{p)/h - eh(k)/;i - [u + Ai/?i])/(eh(k)) 



(85) 



where ep.h(k) and a have been given in Eas. (|63|l and (|74|l . respectively. 

Let us now define the momentum-resolved current, which describes each p process, by 
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<5/(p) ^ / dk|0(q,p)|25(k + p-q) 

X Z(k)5(ef (p)/n - ep(k)/n - + Ai/?i])/(ep(k)) 



+ (a - Z(k))<5(ef (p)/?i ~ eh(k)/?i - [a; + M/?i])/(eh(k)) 



(86) 



We will show that the momentum-resolved current defined in Eq. H8t)|) exhibits a very interesting feature, while we 
may not be able to see any remarkable character by observing the total output coupling current in Ea. l|85|l . The 
physical meaning of Ea. H86() is the weighted current with conserving the total momentum and the total energy. It 
describes two energy-conserving processes. In the first process, the energies of a photon from the external laser field 
and a particle excitation is transferred to a free atom 

ncj + Ai + ep(k) =ef(p). (87) 

In the second process, the energies of photon from the external laser field and a hole excitation is transferred to a free 
atom 

ncj + ^ + eh(k) =ef(p). (88) 

These processes are weighted by the wave function rcnormalization and the Bose distribution function. Since we are 
interested in probing the low-energy excitations in the optical lattice, we evaluate the momentum-resolved current for 
the small momentum p — (0, 0.01, 0.01)7r/a and q = as a function of the energy transfer hu) from the laser field. 
For simplicity, we approximate the low-energy single-particle eigenstate of the free gas by the plane wave, i.e., 

.p(r) ^ ^, .(P) = (89) 



In Figl^l we plot the momentum-resolved current for various temperatures with fixed values of U (Figs|Sta)-(c)), 
and for various values of U with a fixed temperature (Figl^d)). Note the logarithmic scale for the atom current. In 
calculating Ea. (|86|l . we approximate the Wannier function w by the Gaussian wave function [15112^. When U > Uc, 
the energy spectrum of the current has two peaks at any temperature. Obviously, the lower peak is contributed from 
the particle-excitation ep(k) and the upper peak is contributed from the hole excitation eh(k). Thus, the distance 
between the two peaks provides the energy gap Aq. 

From Figsinia)-(c). we find two characteristic features in the temperature dependence of the output coupling 
current. First, the distance between the two peaks increases as the temperature increases, rcfiecting the growth of the 
energy gap in the excitation spectra as shown in Fig|2Ia). Second, the peak corresponding to the particle excitation 
grows as the temperature increases, reflecting the thermally populated particle excitations. One can understand the 
temperature dependence of the particle and hole peaks from the expression for the momentum-resolved current in 
Ea. (|86|l . The particle-excitation peak exists only at finite temperatures since /(ep(k)) ^ as T ^ 0. In contrast, 
the hole-excitation peak exists both at zero and finite temperatures since /(eh(k)) = — 1 — /(|eh(k)|) — > —1 as T — > 0. 
This means that one can only detect the particle excitations and hence the energy gap at finite temperatures. 

In FigEl^d) with the fixed temperature (3zt = 10.0, we find that the particle and hole peaks merge into a single peak 
at U — Uc, where the Mott gap vanishes. As shown in Figl^fa), with fixing U = Uc and increasing the temperature, 
the peak is split into two peaks. This reflects the existence of the energy gap at finite temperatures even at U = [/c as 
shown in Fig|5Jb). This gap, however, is not associated with the Mott insulator phase, as one is in the normal phase 
away from the critical point. 

We now comment on the experimental observability of the particle-hole gap in the output coupling current. In 
most cases the contribution in the current from the particle excitations is much smaller than the contribution from 
the hole excitations. However, the two peaks corresponding the particle and hole excitations are very sharp and do 
not overlap with each other except very close to the critical point. Therefore, the particle peaks should be visible as 
a distinct contribution, although high-precision experiments may be required. 



VI. CONCLUSIONS 



In this paper, we have studied the single-particle excitations in a Bose gas in an optical lattice. We calculated the 
single-particle Green's functions of the Bose-Hubbard model in the Mott insulator phase by using the standard basis 
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FIG. 3: Momentum resolved output coupling current for p = (0, 0.01, 0.01)7r/a and q = 0. Panels (a)-(c) are plots with the 
fixed values of U and various temperatures Pzt = 10.0, 4.0, 3.0, 2.0. Panel (d) is a plot with the fixed temperature with 
U/zt = 5.83, 6.5, 8.0. Here we take A = 852nm and m and as for *^Rb atom, and Vo to adjust the ratio U / zt 



operator formalism whose poles provide the dispersion relations of particle and hole excitations. These dispersion 
relations show that the energy gap between the particle and hole excitations at k = vanishes at the critical point. 
The energy gap increases with the increasing temperature or increasing U/zt. The same behavior of the energy gap 
has also been found in Refs. 0, . 

In order to gain information of the single-particle excitations in the Mott phase, we have calculated the out- 
put coupling current using the linear- response formalism given by Luxat and Griffin 01 . We have shown that the 
momentum- resolved current has two distinct contributions from the particle and hole excitations, and thus one can 
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directly measure the Mott gap. Moreover, this can be used to detect the transition point from the Mott insulator 
phase to the superfluid phase, where the Mott gap vanishes. We note that the contribution from the particle exci- 
tations is strongly temperature dependent, vanishing at T = 0. In most cases the particle-excitation contribution 
is much smaller than the hole-excitation contribution. Nevertheless we expect the particle-excitation peak to be 
experimentally observable since it is distinct from the hole peak, as shown in Fig|31 

In the present work we have restricted our calculations to the Mott insulator phase. It is next important to study 
how the output current changes as we enter the superfluid phase. In the superfluid phase, we expect to have a direct 
tunneling current from the condensate component, as well as processes involving quasi-particle excitations jHSH- 
Extending the present study to the superfluid phase will be give n in a future work. 

Finally, we mention the recent experiment by Stoferle et al. p30l | investigating excitations in a trapped Bose gas in an 
optical lattice using the Bragg spectroscopy. This type of experiment deals with the density-fluctuation spectrum P, 
l3lll32l|. In contrast, our present paper probes the single-particle Green's functions in a Bose gas in an optical lattice, 
which is more directly related to the single particle and hole excitations that contain important information about 
the superfluid-Mott insulator transition. 
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APPENDIX A: DERIVATION OF THE LANDAU FREE ENERGY BY USING THE INVERSION 

METHOD 

This Appendix gives a detailed derivation of the Landau free energy for the Bose-Hubbard model used in Sec. ^1 
Our aim is to derive a free energy as a function of an order parameter describing the superfluid phase, defined by 
= (bi), by perturbative expansion in the hopping term. Following Ref. (see also Refs. |33,I3), we introduce 
an expansion parameter g and write the Hamiltonian as 

H{g)=Ho{9) + gH,, (Al) 

where 

i i 

-J2^vn9)b^ + mi9)b% (A2) 

i 

Til = ~tJ2 blbj. (A3) 

<i,j> 

In Ea. (|A2p . we have introduced the symmetry-breaking external fields 77 and 77*, which controls the order parameter 
^'i. We regard 77 and 77* as functions of the coupling constant g, imposing the condition that the relation 

should hold independent of g. 

The thermodynamic potential associated with the Hamiltonian Ea. (|Al|) is given by 

ri[77 , 77*] = -ilnTr e-^^(9). (A5) 

P 

By performing the Legendre transformation, the Helmholtz free energy as a function of the order parameters 4" and 
4** is given as 

^[^,^*]^n[^^,lf] + Y,[vK9)'^^ + v^{9m]■ (A6) 

i 

Differentiating the above equation in regard to we find 

^^^^^V^[9■,^^,^*]. (A7) 
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The required free energy is obtained by setting g = 1 and turning off the artificial external fields, i.e., r]i{g = 1) = 
at the end of the calculation. Then, the equation that determines is given by 



dm 



= 0. 



(A8) 



We now expand the thermodynamic potential in powers of g by using 

g-/30[g;r,,r,*] ^ rj,^. g-/3W(s) 



= E(-f) A / ^^i-- - / '^^"Tr {e-^'^oWT. [Hi(ri) • • • 7^i(r„)]} 



= Tr e"'^'^"(") 



dri Tr 



H(ri) 



dri / rfT2Tr 



-/3Wo(0) 



H(ri)Hi(r2) 



(A9) 



where Hiir) = e^'^'^^'^/'^Tiie ''^o(o)t/?i^ g^^-^^-j ^-j-^^ operator T^ is the imaginary time ordering operator. We expand 
^[9;V,V*] as 



1][.9; ^7, r;*] = f!'") [ry, ,7*] + gf]^) [,/, rj*] + g^f]^^) + . . . . 

In our formalism, 77 is also expanded through the relation ljA4|l as 

r][g; ^, **] = + 5?/^^ ^'*] + g^^/z) + 

By the Legendre transformation, the Helmholtz free energy is expanded as 

r[^, m] = r(°) m] + gv'-^^ m] + y^r^^) [^-, ^f*] + . . . , 

The zeroth-order term F is given by 

i 

where rj^^^ is related to ^ through (jA4|l with ^ = 0. The next two terms are given by 



(AlO) 
(All) 
(A12) 

(A13) 



E 

i 

J7(l)[^(0)^^*(0)] 
J7(2)[^(0)^^*(0)] 



I 9^ 



(0) 



(0) 



+ H3 



dr]i 



(AM) 



(A15) 



Following the procedure described in Ref. oii^ <^an successively calculate Ty*-"-*, and hence systematically expand 
r in powers of the coupling constant g. 

We now explicitly derive the free energy as a function of the order parameter to first order in the coupling constant 
g. The zeroth-order term of the thermodynamic potential is given by 

O(0)[^(0),^*(0)] ^ -^InTr e-^«(°) 



InTr 



(A16) 
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where we have defined the local Hamiltonian as 

^ ^blblhh - ^^^blh - {vf'>b. + vl'^b]) . (A17) 

As noted above, rif^^ is a function of ^'i through the relation l|A4p . but it is in general difficult to obtain an explicit 
functional form of rj^^^'^i, 5"*]. In the spirit of the Landau free energy, we assume that the order parameter is small 
and thus we expand Qi in powers of rjf^'' and jy*'"' as follows: 



+ {{h.) + ^M)-{m]\4'^\\ (A18) 



where we have introduced / and h as 



, _ n+l n CAIQ^ 



h 



(n + l)(n + 2) (n+l)2 



{Un - ^l,)^Ui2n + 1) - 2^l,] {Un - /i,)^ 

n{n + 1) 



[U{n - 1) - (f/n - ^HnU{n - 1) - /i- 

?i(rt — 1) 

nin + 1) 



{Un- iMi)[U{n-l)- 
and defined the local thermal average by 



with 



(A20) 



Zo^^e-^^", (A21) 

n 

(••■)^^E(---)^"'''"' (A22) 



En = ■^"'('^ ~ 1) ^ (A23) 

Differentiating fi^"' [t?^"' , ?7*("'] with respect to 77*, we obtain = [77'"' , 77*^°)] 

^ _ 9r)(0)[,7(0),r;*(")] 

= (/)7yf + {2(.9(n))+/3[(/2)-(/)2]},7f P 

= *f [77(°),77*W]. (A24) 
The external fields are expressed by the order parameters by inverting Eq. 1A24(I as 

ryf'[*] = (A25) 
Next, we calculate the 1st order term of the Landau free energy. From Eq. I)A14|I . we find 

r(i) **] = (Wi),=o = -i E (^^^.%=o, (A26) 

<i-J> 

where the expectation ( )g=o is evaluated with respect to 7i(0). We find that 

r(i) = -t ^ (A27) 

<i,i> 
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Then 77^^^[*, ^*] is given as 



Up to first order in the couphng constant we have 



^ (0) I (1) 
.(0) 



(A28) 



(A29) 



Turning off the external field, i.e., setting 77 = 0, we find ry- = ~iX]<j>'^j- -^y substituting this expression to 
Ea. (jA25|l . we can obtain a closed equation that determines We note that our first order calculation is equivalent 
to the mean- field theory in Ref. [HIS, since 'qf'^ plays a role of the internal mean field obtained from the decoupling 



h\h,^{b\)h, + {h,)b\ = %b,+^.,b]. 



(A30) 



The advantage of the present formalism is that one can systematically include the higher-order terms in the internal 
field r?f^ 

By using Ea. HA13|) . (|A14|I and (|A29|I . we finally obtain the expression for the Helmholtz free energy to 1st order in 
the hopping parameter t. For simplicity we present the result for the homogeneous case. The free energy per lattice 
site is given by 



\nZo + zt{l~zt{f))m^ 



(h) + I: m {f? 



Comparing Ea. lpOT)) with Eq.Q, we find T^") = - lnZo//3 with defined by llMT|l and 

A = zt{l-zt{f)), 

B = (zt)4{(/.)+/3[(.n-(/)2]}, 



with / and h defined by Eas. (|A19p and (|A20p . At T = 0, these coefficients reduce to 

n + 1 n 



A = zty.- zt 
B = {ztf 



Un — ii C/(n — 1) — /i 
{n + l){n + 2) 



(A31) 



(A32) 
(A33) 



(A34) 



(n+1)^ 

([/n- ^)2[[/(2n+ 1) - 2/i] {Un-^y 
n? n{n + 1) 



n{n — 1) 
[U{n - 1) - fj]^[Ui2n - 3) - 2^] 

n{n + 1) 
(/7n-/i)[C/(n- 1) 



(A35) 



These results for T = agree with Ref. [l^. The phase boundary between the Mott insulator (or normal) and 
superfluid is determined by A = 0, and thus zt{f) = 1. 



APPENDIX B: GREEN'S FUNCTIONS 



The four components appearing in the expression for Gi(k, w) in Ea. (|57|l are given by 

G(„)(,i+i),(Ti+i)(„)(k, w) 



hPn,n+l (jT-i^ - Vn + K-1 " -Pn,n-l7'(ti)(„_i), („_!)(„)) 



{huj - ep(k)) {fiuj - e/,(k)) 



(Bl) 
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G'(„)(ri+l),(ri)(ji-l) (k, Uj) 
pk 

-(n+l)(n),(n-l)(n) 



G(„_i)(„)^(„+i)(„) (k, uj) 

^ ^^n,n+l-Pn.n-lT'(^)(„_l) („)(„+!) 

- ep(k)) (ritj - e/i(k)) 
G'(„_i)(,i),(„)(„-i) (k, u) 



(ri+l)(n),(n)(n+l) 



{huj - epCk)) {huj - ehCk)) 
while the four components appearing in the expression for G2(k, cj) in Ea. (|58|l are given by 

G(„+i)(,i),(Ti)(„+i)(k,a;) 



(atj + ep(k))(?iw + eft(k)) 
G(„)(„_i), („)(„+!) (k, uj) 

^ ^Pn+l,nPn-l,nT^n-l){n),(n+l){n) 
G'(n+l)(n),(n-l)(n) (k, Uj) 



T^Pn , n - 1 Pn , n + 1 Pf^ 



(n)(n+l),(n)(n-l) 



G(„)(,i_i),(„_i)(,i) (k, uj) 

^Pn,n-l ( — + Vn+1 — Pn,n+lTf^ 



(n)(n+l),(n+l)(n) 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



{huj + ep(k)) {huj + eh(k)) 

where, Da — Djs = Pa,p and Va^a = K,, and the expressions for T and are given in Eas. l|29|l and (|30|l . respectively. 
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